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Abstract . 

Some ideas for using hodograph theory, mapping techniques and method of characteristics to formu- 
late typical aerodynamic design boundary value problems are developed. Inverse method of charac- 
teristics is shown to be a fast tool for design of transonic flow elements as well as supersonic flows 
with given shock waves. 



Introduction 

This paper is intended to illustrate a revitalization of classical tools of theoretical aerodynamics for 
use on modem graphic workstation computers presently available to the design engineer: 

The theoretical methods have their origins in the time before large scale numerical computing be- 
came the standard approach for analyzing aerodynamic performance, about two decades ago these 
tools were already operational for practical aerodynamic tasks. Transonic flows then posed challeng- 
ing problems and analytical mathematical modelling was used to gain insight into various theoretical 
and applied questions resulting from nonlinear model equations. The Hodograph Method gave an- 
swers to many of such type problems. Similarly, the “Method of Characteristics” yielded practical re- 
sults. Combination of both methods, more recently, has permitted inverse - or at least indirect - 
formulation and solution of aerodynamic design problems. 

Later, numerical methods became more important because of their general applicability but frequent- 
ly they give only poor insight into a mathematical model underlying a described phenomenon. The 
previous hodograph approach was somewhat complicated because of mapping procedures, but results 
still serve as test cases for numerical methods. 

Nowadays, while most of the successful analytical methods are used only for educational purposes, 
we witness another type of tool emerging from developments in computer technology: Graphic 
Workstations and even PC’s provide powerful computation, illustration and documentation of results 
to the model equations for fast aero analysis and design. Interactive methods are being developed to 
provide a strong coupling of computer power and speed with the design engineer’s experience and 
strategies to obtain his design goals: For many applications we seem to have the knowledge base and 
computer hardware to develop a variety of what may be called Aerodynamic Expert Systems . 
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In this situation we may want to recall classical “pre-CFD” methods because, if implemented to those 
ast graphical deslctop computers, they might be modernized and improved to easily give fast first 
steps for aerodynamic design and optimization, and last not least to serve as educational tools. 

In the present paper we illustrate the idea of combining fast classical aero methods with most recent 

— ,T d T teChn ° l0gy by using hod °graph formulations and characteristics to obtain 
ome well known and some new plane and axisymmetric transonic and supersonic flow elements 
The fast computation and powerful graphic evaluation of results invite experimenting with conceptu- 
al extensions. Here the hodograph method is extended to axisymmetric flows and a method of char- 
acteristics for axisymmetric rotational flows will be presented and proposed for use of designing 
more general three-dimensional flows. ® 


Hodograph-based methods for transonic flows 

“ Wing 7 view of .a" > extended hodograph method is focused on transonic applications. For 

c!^ n ? f 10na 1Sen fl tr0piC ^ ow thls a PP r °ach is wellknown in the literature, here the illustration is 
earned out for plane flow and extended to small perturbation axisymmetric transonic flow - an option 

foSsy^ertc fbw USC C main PUrP ° Se ° f thC h0d0graph ’ Iinea nty, obviously cannot be obtained 
Potential flow models 

Isenmipic flow assumptions result in a system of PDFs for potential <D and streamfunction V with D 
v^riableT 135 ^ dimenS1 ° nleSS density ajld velocity, and flow angle d and velocity Q independent 

W 2 - !) '*'* 


d> = 

* D Q 


( 1 ) 


Rheograph transformation: Beltrami equations 


Here we use a modification of the hodograph variables: System (1) is transformed 
Prandtl - Meyer angle 


by using the 


v = J V|a/ 2 -i|^ 

(Q= l) U 


( 2 ) 


as one independent variable instead of the velocity Q. The basic PDEs become 
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*v = JK*t 


<f> = 


KY. 


( 3 ) 


with j = -1 representing subsonic flow (where v < 0) and j = 1 for supersonic flow (where v > 0). The 
coefficient 


K = 



(4) 


is, for isoenergetic flow and with (2), a function of v only. So far any hodograph problem formulation 
is just a matter of stretching from a (Q, ■&) - plane to a (v, fl) - plane. 

The technique stressed here involves further elliptic or hyperbolic mapping of the variables (v, i3) 


v 

S 



= J*S 


which results in a transformation of (3) to become 


* = J™, 


<J>, = KV 

t S 


(5) 


( 6 ) 


Both (3) and (6), and also the Cauchy-Riemann or wave equations (5) are more generally named Bel- 
trami equations. With K now a function of s and t this system of PDEs is linear. The first author has 
made extensive use of it for transonic airfoil design [1, 2]. 

For the purpose of illustrating a generalization to axisymmetric flow, we use here the small perturba- 
tion version of (5) and (6), with a notation 


v ~ U 


t - VY~ Pl 


<t>~X 


(7) 


__ yl + Pt 


where for plane 2D flow pi = 0. Replacing the coefficient K by its leading term near sonic conditions, 
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jK~^ (g) 

with p 2 = 1/3 results in the near sonic version of the systems (5) and (6), the compatibility relations 



K = ^U s 

and characteristic equations 

x g = u Pi y 1 


(9) 


Now we see an elegant symmetry of these two coupled pairs of equations, each modeling a general- 
ized axially symmetric potential, [3]. We can distinguish between various types of flow, depending 
on the parameters j, Pjand p 2 . Linear subsonic (j = -1) or supersonic (j = 1) flow is described by p? = 
0, while transonic flow requires p 2 = 1/3, with j both -1 and +1 for mixed type flow. With Pl = 0 or 1 
we have plane 2D or axisymmetric flow, respectively. Mapping in various aero or fluid dynamics 
case studies can so be reduced to one generalized system of basic equations [4, 5], Any one of p, or 
P 2 being equal to zero yields linear equations, but for near-sonic axisymmetric flow a weak nonlin- 
earity persists, which seems to be the reason why this formulation has not been used for aerodynamic 
problems, except in the one work by Hassan [6], 



Fig. 1. Rheograph or Characteristics plane: Elliptic (shaded) and hyperbolic (cross-hatched) 
domain for mixed type model equations 


Self-similar solutions 

We see the relation of Beltrami equations to conformal and characteristic mapping: singular solutions 
in classical hydromechanics have helped to understand many aerodynamic phenomena, so we wish to 
use the system for axisymmetric near-sonic flow also for solving some of its typical features. 

Figure 1 illustrates the working plane (s, t): neither physical plane (X, Y) nor hodograph plane (U, 
V), it is suited for a definition of boundary and initial value problems which require a parametric for- 
mulation. In transonic flows, the mixed elliptic/hyperbolic type subdomains require contact along the 
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mapped sonic line, here suitably fixed at s = 0. Earlier applications and illustrations [5] explain the 
use of the name “Rheograph” and “Characteristics Plane”. 

Some classical and many new phenomena may be modeled from the general harmonic set of self- 
similar solutions in polar coordinates 


U = r n ■ h (<p) 


V = r n+Pl b k(y) 
X = r b+Pl ' n -f(<p) 


( 11 ) 


Y = r b ■ g (tp) 

which require only solving a set of four coupled ODEs for the generalized harmonic functions h, k, f 
and g, with two free parameters n and b. 



Fig. 2: Quasi-harmonic functions for far-field singularity U,V(X,Y) modeling flow past a body 
of revolution in sonic free-stream Mach = 1 


Example: Guderley’s far-field. singularity of an axisymmetric body in sonic flow 

M. Klein [7] has investigated these coupled potential flow problems calculating some plane and axi- 
symmetric cases with different exponents n, b. This was done prior to using some gained knowledge 
for setting up more general boundary/initial value problems for numerical solution of (9) and (10) 
with a Poisson solver and the method of characteristics on a graphic workstation. One case studied in 
detail is the solution for simulating the flow past a body of revolution in sonic free-stream. This is a 
classical transonic problem first solved by Guderley 1954 [8] and elegantly confirmed by Muller & 
Matschat 1964 [9]. Their work suggests use of (1 1) with a ratio of the exponents b/n = - 7/9. In fact, it 
is just this ratio which yields a physically reasonable solution. 
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Figure 2 illustrates the result: Graphic CFD postprocessing software is used to show the U and V dis- 
tribution in the physical meridional plane (X, Y). This example may be used to illustrate the use of 
working in the Rheograph plane, to understand flow details with nonlinear model equations better 
and to have more freedom to suitably model boundary values. In transonic flow, meaningful solu- 
tions frequently can only be obtained by formulating boundary conditions in an indirect, inverse way 
- this is the basic reason why some practical design problems are easier solved in inverse mode. 

For the following transonic 2D example (pj = 0) we return to the systems (5) and (6), the Rheograph 
equivalent of the 2D full potential equation. 

Example: 2D transonic nozzle exit 


Equations (5) and (6) for supersonic flow j = 1 transform into compatibility relations 

dd\ 


dv 


= 1 


rj = const 


and characteristic equations 


d-d 



= const 


= -1 



const 


_1 

K 


dV 


4 


= const 


K 


( 12 ) 


(13) 


which are the basis for a rapid linear method of characteristics. Implemented on a graphic worksta- 
tion, solutions may be obtained and visualized extremely fast, we use the method to set up a knowl- 
edge base for interactive transonic design expert systems with advanced graphic pre- and 
postprocessing. Flexible geometry input for boundary conditions was used by Gentner [10] to define 
a 2D sonic throat and the downstream accelerated exit flow. With initial data for Mach number flow 
angle and physical coordinates along the t-axis (Fig. 1) a first calculation determines the solution of 
(12) and with K(v(^,T])) available, the second step is the solution of (13). 

Figure3 once more stresses the difference between hodograph and Rheograph or Characteristics 
plane: The flow structure may map into a multivalued hodograph, while the Rheograph may be con- 
trolled to show a single-valued characteristics grid. 

The result with a non-symmetrical exit contour designed by prescribing velocity distribution along 
the nozzle axis is depicted in Figure 4. The idea here was the combination of (known and well-devel- 
oped) potential flow modeling with mapping transformations based on hodograph theory (also 
known but considered complicated), and the use of powerful workstations (helping with rapid com- 
putation and graphic visualization). The above isentropic model equations are either linear or weakly 
nonlinear. In the following, design problems involving non-isentropic flow will also be solved by the 
method of characteristics. 



Fig. 4: Laval nozzle exit designed from sonic line Cauchy data and velocity distribution (Mach, 
flow angle) along a curved axis. Color graphics illustrate stream function / contour . 






289 

Tnirc »ru*mot*ono: CoMerwx:* <y- wvtrw Cooc*pts ooc Oottmizouon in Erujtn^rmg Sc*«nc*5 

f ICOEI5- H! , Ec , :D f .C.S.Ouii^ r ovic^ v»osnmg*o*" C.C..C ) c*ope r 22 - 21 . 1 99 ' 


Numerical methods of characteristics have been developed for plane and axisymmetric , for isentrop- 
ic and rotational flows. Prescribing arbitrary shock waves results in rotational flow because of shock 
curvature. Cauchy initial data for flow field computation therefore require coordinates, velocity com- 
ponents and entropy distribution along the prescribed shock geometry. For Cross-Marching (Fig. 6b), 
the iterative calculation of entropy convection along the streamlines requires an extrapolation of data 
BD— >C, while for the usual downstream marching (Fig. 6a) an interpolation of data A— >D<— B is 
needed. 

The following two examples were obtained with a new numerical Cross-Marching method of charac- 
teristics for axisymmetric isentropic or rotational flow by the second author [ 12]. A flexible input ge- 
ometry generator and workstation implementation lay ground for further extensions and use for 
aerodynamic design tasks. 

Example: Segment of a conical flow field 

As a first example for the new method of characteristics a part of the flow field past a circular cone is 
computed. Input data are the upstream Mach number, a set of coordinates of and post-shock condi- 
tions behind the given conical shock wave with given angle. Fig. 7 illustrates the characteristic grid, a 
selected integrated streamline and reveals a limit line singularity along a ray through the cone vertex, 
well within the solid cone which is compatible with the shock cone and Mach number. The case is 
well-suited for checking the accuracy in comparison with the solution of the Taylor- Maccoll ODE; 
graphic visualization of the Mach number and flow angle distribution must show constant values 
along rays through the cone vertex, though its location is not part of the input data. 



Fig. 7: Ideal gas (y = 1.4) flow past a circular cone. Given Mach = 2, shock angle = 45° (com- 
patible with a solid cone of 27.32°). Choice of shock segment size relative to axial distance: 
Characteristics grid (a) with or (b) without limit cone of -16.3°. Every third grid line shown. 
Color graphics (iso - Mach) for flow field conicity check (c). 

Computation time on a Sun Sparc Station: 6 seconds. 
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Example: Segment of a flow field downstream of a curved shock 

A slight variation of the input shock segment geometry brings rotation to the flow field downstream 
of the shock. Color graphic visualization of the velocity and pressure distribution shows a strong de- 
viation from conical structure, Fig. 8. 





Fig. 8: Flow with a curved shock wave. Given Mach = 2, shock angle varies from 50° to 45°. 
Characteristics grid, surface streamline (a). Color graphics for iso-Mach (b) and iso-flow angle (c) 


Design of three-dimensional flow fields 


The exploitation of plane and axisymmetric inviscid flow fields for the definition of flow patterns 
generated by three-dimensional bodies in supersonic flow has been used since about three decades 
when Nonweiler [13] created the first “waveriders”. In recent years renewed interest originated in 
such configurations for generic lifting aerospace transport vehicles and supersonic inlet shapes [14]. 
The first author recently contributed an idea to this research which is aimed in generalizations of 
waverider shape definition by applying conical flow solutions with constant shock strength but axial 
distance of the shock segments varying along span [1 1], The idea is based on the assumption of a lo- 
cal axisymmetry in every 3D flow, which is well defined in an”osculating plane” if the shock wave is 
known. Fig. 9. The method to generate three-dimensional configurations requires little more effort 
than evaluating one Taylor-Maccoll conical flow solution. Only the Mach number, shock angle, lead- 
ing edge and shock profile in the exit plane need to be prescribed. Based on this method a very rapid 
interactive design code was developed by Center et al [15]. Numerical analysis with an Euler code 
shows a striking agreement of the numerically captured shock location with the design solution 
which makes it worthwhile to further develop such techniques. 

The method of characteristics with Cross-Marching developed here may be used for such further de- 
velopment, with the possibility of exploiting also rotational flow fields to find 3D body contours 
This method is then equivalent to solve the 3D Euler equations in an inverse design mode, rapidly 
earned out in an interactive fashion on the workstation. 
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Fig. 9: Local conicity in osculating plane of a flow behind given rule surface shock wave: Design 
of generalized super/hypersonic waveriders. 


Conclusion 

We have tried to illustrate some ideas to extend classical theoretical methods for the aerodynamics of 
inviscid, compressible flows. The purpose is an implementation of these tools to develop software on 
fast modem workstation computers which enables the design aerodynamicist to perform rapid early 
stage design studies with aerodynamic expert systems, but also to develop these techniques toward 
convincing educational programs for students. Transonic and supersonic aerodynamics require in- 
verse problem formulations if flow properties should be used optimally for design goals. This was 
shown using the method of characteristics for inverse applications. 
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